Evidence for a Massive Protocluster in S255N 



C.J. Cyganowski 1 , CL. Brogan 2 , T.R. Hunter 2 
ccyganowOastro . wise . edu 
cbrogan@nrao . edu 
thunterOnrao . edu 

ABSTRACT 

S255N is a luminous far-infrared source that contains many indications of 
active star formation but lacks a prominent near-infrared stellar cluster. We 
present mid-infrared through radio observations aimed at exploring the evolution- 
ary state of this region. Our observations include 1.3 mm continuum and spectral 
line data from the Submillimeter Array, Very Large Array 3.6 cm continuum and 
1.3 cm water maser data, and multicolor IRAC images from the Spitzer Space 
Telescope. The cometary morphology of the previously-known UCHn region 
G192. 584-0. 041 is clearly revealed in our sensitive, multi-configuration 3.6 cm 
images. The 1.3 mm continuum emission has been resolved into three compact 
cores, all of which are dominated by dust emission and have radii < 7000 AU. 
The mass estimates for these cores range from 6 to 35 M . The centroid of the 
brightest dust core (SMA1) is offset by 1.1" (2800 AU) from the peak of the 
cometary UCHn region and exhibits the strongest HC3N, CN, and DCN line 
emission in the region. SMAl also exhibits compact CH 3 OH, SiO, and H 2 CO 
emission and likely contains a young hot core. We find spatial and kinematic 
evidence that SMAl may contain further multiplicity, with one of the compo- 
nents coincident with a newly-detected H 2 maser. There are no mid-infrared 
point source counterparts to any of the dust cores, further suggesting an early 
evolutionary phase for these objects. The dominant mid-infrared emission is a 
diffuse, broadband component that traces the surface of the cometary UCHn re- 
gion but is obscured by foreground material on its southern edge. An additional 
4.5 /iin linear feature emanating to the northeast of SMAl is aligned with a 
cluster of methanol masers and likely traces a outflow from a protostar within 
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SMA1. Our observations provide direct evidence that S255N is forming a cluster 
of intermediate to high-mass stars. 

Subject headings: stars: formation — infrared: stars — ISM: individual (S255N) 
- ISM: individual (G192.60-MM1) — ISM: individual (G192. 584-0.041) — tech- 
niques: interferometric — submillimeter 



1. Introduction 



The formation and e arly evolution of ste llar clusters occurs in deeply embedded regions 
of giant molecular cl ouds (ILada fc Lad a 2003). While much has b een learned from recent sur- 
veys in the infrared ( Gutermuth et al. 20051 : Muench et al. 2002 ). the earliest stages of clus- 
ter formation will (at least in many cases) be hidden from all but the longest IR/millimeter 
wavelengths. Due to the small size scales of youn g clusters-multiplicity of protosta rs has 
been observed on scales of a few thousand AU (e.g. iMegeath. Wilson, fc Corbinll2005l )-high 
angular resolution is necessary to resolve individual objects, particularly in massive star 
forming regions which lie at relatively large distances (> 1 kpc). These considerations point 
to millimeter- wavelength interferometric observations of thermal dust continuum emission as 
an effective means of searching for clusters of young protostars, as the mm emission remains 
optically thin up to high column densities (]Vh~10 25 cm -2 ). 

Located ~ 1' north of the luminous infrared cluster S255IR, S255N is a promising tar- 
get in the search for young protoclusters. As illustrated in Figure (H S255N is located in 
a complicated region of past and ongoing massive star formation. S255N and S255IR (sat- 
urated in this mid-IR Spitzer IRAC image) lie between two large Hn regions, S255 and 
S257. Large-scale 12 CO and HCN observations show that S255IR and S255N occupy op- 
posite ends of a molecular ridge between the Hn regions (IHever et al.lll989l). The total 
luminosity of S255N (~ 1 x 10 5 L ) is about twice that of S255IR jMinier et alil2005h . and 
single-dish continuum and spectral line observations at submillimeter and millimeter wave- 
lengths have established the presence of large column den si ties of dust and gas toward both 
regions fe.g. iRichardson et al.l Il985l : iMezger et al.l 1 19881 : IZinchenko. Henning. fc Schreyer 



19971 : iMinier et al.ll2005l ). Observations at infrared and radio wavelengths, however, sug- 
gest that S255N is the younger of the two regions. For example, S255IR is bright (> 
70 Jy) in all infrared bands of the Midcourse Space Ex periment (MSX) and contains a 
well-developed n ear-IR cluster of early-B type stars fS2 55-2lHoward. Pipher fc Forrestll 19971 : 
2001|) , a cluster of co mpact Hn regions JSnell fc Ballvl Il986l ). and a wealth 



Itoh et al 



of complex H2 emission features (IMiralles et al.l 119971 ). In contrast, S255N (also called 



S255 FIR1 and G192.60-MM1) is undetected in MSX images at wavelengths shorter than 
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21 yum (ICrowther fc Contil l2003|j. and contains only a single cometary UC Hll region, 
G192.584-0.041 (e.g. lKurtz. Churchwell. fc WoodlEoSh . 



Additional evidence for protostellar activity in S255N exists in the form of outflow 
tracers. For example, two small knots of H 9 emission bracket the UC Hll region and may 
be tracing an outfl ow (IMiralles et al.lll997l ). In beamsizes o f ~ 50", redshifte d CO emission 
( iHeyer et al.l Il989l ) and highly-blueshifted OH absorption (jRuiz et al.l Il992l ) are also seen 
toward the UC Hll region. Finally, 44 GHz (7o-6i) A + Class I methanol masers form a 
linear feature extending northeast of the UC Hll region (jKurtz. Hofner. &: Alvarez! I2004J ) ; 
masers of this type have been observed in association with m olecular outflows in other objects 
(Ipiambeck & Mentenl [l99ol : iKurtz. Hofner. fc Alvarez! l2004h . 



At this stage, further understanding of the state of star formation in S255N requires 
resolving the mm dust continuum emission in order to search for additional compact sources 
that may be present in the vicinity of the UC H II region. The high angular resolution now 
available with the Submillimeter Array (SMA)0 makes this goal possible, while the wide 
bandwidth allows simultaneous observation of many spectral lines, which are sensitive to a 
range of gas conditions across the region. We describe our observations in section 2, present 
our results in section 3, and discuss our interpretations in section 4. 

A ra nge of distances to S255 and S255N have been used in the literature. At the 



extremes, 



Georgelin. Georglein. fc Rouxl (119731 ) find an optical distance of 3.2 kpc to each 



of the adjacent Hll regions S254 and S257 (located west of S255N), while iHunter fc Massey 



(119901 ) find a distance of 1.1 kpc to the optical Hll region S255 (located east of S255N, see 
Figure [[]) based on spectroscopy of its exciting star. Using a LSR velocity of +9 km s _1 
(typical of the centroid velocities of our observed lines), and a Galactic center distance 



Brand & Blitz! Jl993h. In this 


paper 


Moffat et al. 


1979|; 


Minier et al. 


2005) 



lr The Submillimeter Array is a joint project between the Smithsonian Astrophysical Observatory and the 
Academia Sinica Institute of Astronomy and Astrophysics and is funded by the Smithsonian Institution and 
the Academia Sinica. 



-4- 



2. Observations 

2.1. Submillimeter Array (SMA) 

Our SMA observations toward S255N were obtained on December 06, 2003 in the com- 
pact configuration. Six antennas were operational and the projected baseline lengths ranged 
from 12 to 62 kA, resulting in a synthesized beam of 477 x 2"4 (P.A.=— 45.65°). In this 
configuration, the interferometer is not sensitive to smooth structures larger than 17". The 
phase center was 06 h 12 m 53 s .56, 18°00'25".0 (J2000), and the double-sideband SIS receivers 
were tuned to an LO frequency of 222.78 GHz. With a bandwidth of ~2 GHz, the cor- 
relator covered 216.796-218.764 GHz in lower sideband (LSB) and 226.796-228.764 GHz in 
upper sideband (USB). The data were resampled to provide uniform spectral resolution of 
1.12 km s . The zenith opacity as measured by the 225 GHz tipping radiometer at the 
Caltech Submillimeter Observatory (CSO) was 0.18 at the beginning of the observations 
and fell as low as 0.16. Data recorded late in the observing session, when the opacity had 
climbed to 0.3, were not used. The typical system temperature at source transit was 200 K. 
The primary beamsize of the 6-meter SMA antennas at this frequency is 56". 

Initial calibration of the data was performed in Miriad. The gain calibrators were 
J0739+016 and J0423-013. J0423-013 was also used for passband calibration. Subsequent 
processing was carried out in AIPS. The continuum emission was measured using line-free 
channels and removed in the u — v plane. The resulting continuum-only data were then self- 
calibrated; these complex gain solutions were also applied to the continuum-subtracted line 
data. The absolute position uncertainty is estimated to be 0."3 and the amplitude calibration 
is accurate to 20%. For maximum sensitivity, both the continuum and line data were imaged 
with natural weighting. The rms sensitivity of the continuum image is 4 mJy beam -1 (8 
mK) and the rms sensitivity of the line data is 89 mJy beam -1 (170 mK). 



2.2. Very Large Array (VLA) 

Archival 3.6 cm data from the NRAo| Very Large Array (VLA) were calibrated and 
imaged in AIPS. The observation date was 2003 June 15 and the total time on source 
was ~70 minutes (project code AH819). The VLA was in A-configuration, in which the 
interferometer is not sensitive to smooth structures larger than ~7". The bandwidth was 
50 MHz in two IFs. The flux calibrator was 3C147, and the gain calibrator was J0613+131. 



2 The National Radio Astronomy Observatory is a facility of the National Science Foundation operated 
under agreement by the Associated Universities, Inc. 
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The synthesized beam of the 3.6 cm continuum image is 0V27 x (T/23 (P.A. =77.51°) and the 
rms sensitivity is 18 /xJy beam -1 . Archival 3.6 cm data from the B configuration (in which 
the interferometer is not sensitive to smooth structures larger than ~20"), observed on 1990 
August 27 (project code AM301), and from the the CnB configuration, observed on 2005 
June 21 (project code TSTCJC), were also reduced and combined with the A-configuration 
data. The resulting image, made with a UV taper at 500 kA, has a synthesized beam of 
1"03 x 0"84 (P.A.=-81.42°) and a rms sensitivity of 34 fiJy beam -1 . 

Archival 1.3 cm data from the VLA A-configuration (project code AC299) were analyzed 
for water maser emission at 22.235 GHz. The observation date was 1991 August 1, total 
on-source integration time was ~10 minutes. The phase center of this observation was 
toward the IR cluster ~ 1' to the south of S255N, hence a correction for the primary beam 
attenuation has been applied to the data. The bandpass calibrator was 3C84 and the flux 
calibration was derived assuming a flux density of 2.16 Jy for J0528+134. The synthesized 
beam is 0'/18 x 0'/16 (P.A. =—56.55°), the spectral line channel spacing is 0.33 km s -1 , and 
the rms noise is 0.13 Jy beam -1 . 



2.3. Spitzer Space Telescope 



Mid-infrared images of S255 were obtained with the IRAC camera (IFazio et al.ll2004l ) 
on the Spitzer Space Telescope as part of Guaranteed Time Observations program 201 (P.I. 
G. Fazio) on 12 March 2004. Integrations of 0.4 s and 10.4 s were taken in the high dynamic 
range mode; S255N is not saturated in the longer exposures, and only the 10.4 s exposures 
are discussed in this paper. Four 10.4 s exposures covered S255N, for a total integration time 
on S255N of 41.6 s. Mosaiced post-BCD 3.6, 4.5, 5.8, and 8.0 /im images, calibrated and 
processed using pipeline version S13.2.0, were downloaded from the Spitzer data archive. 



2.4. Caltech Submillimeter Observatory 



Our submillimeter continuum observations were obtained at the CSO using the Submil- 
limeter High Angular Resolution Camera (SH ARC), a 3 He-cooled monolithic sil i con bolome- 



ter ar ray of 24 pixels in a linear arrangement (IHunter. Benford & Serabynlll996l ; I Wang et al. 
19961 ). For a typical dust source with a submillimeter spectral index of ~ 4, the effective 



frequency of the broadband 350 /im filter is 852 GHz and the bandwidth is 103 GHz. An 
on-the-fly (OTF) map of S255 was obtained on 21 December 1995 by scanning the array 
through the source in azimuth while the secondary mirror was chopping at a rate of 4.1 
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Hz and a throw of 88". Successive scans were made after stepping the array in elevation 
by increments of 5". Airmass corrections were applied to each scan using the opacity de- 
rived from frequent scans of Saturn duri ng the night. The map data we re restored with a 
N0D2 dual beam restoration algorithm (jEmerson. Klein &: Haslamlll979l ) and transformed 
into equatorial coordinates. The resulting image was smoothed with a Gaussian to produce 
an effective half-power beamsize of 15". 



3. Results 



3.1. Continuum emission 



Our 1.3 mm SMA continuum data resolve three distinct sources withi n the previously- 



observed submm/mm clump of S255N (aka S255 FIR1 and G192.60-MM1: Ijaffe et al.lll984 



Mezger et aDll988l : IMinier et al.ll2005h . Figure [2] shows the SMA 1.3 mm continuum image 
with CSO 350 fim contours superposed. As illustrated in Figure [2J the strongest SMA 
1.3 mm emission peak coincides with the CSO 350 /zm peak (resolution 15"). The CSO 350 
/im integrated flux density is 575±20 Jy, consistent to wit hin 10% of the val ue predicted 
by the dust spectral energy distribution (SED) models of IMinier et al.l (120051 ). The three 
mm sources resolved with the SMA are designated SMA1, SMA2, and SMA3 in order of 
descending peak intensity. The observed properties of each source (peak intensity, brightness 
temperature, integrated flux density, and source size) are listed in Table [TJ and the sources 
are labeled in Figure [2J The integrated flux densities and source sizes listed in Table [TJ were 
determined by fitting a single Gaussian component to each SMA source. SMA1 was not well 
fit by a single Gaussian, indicating that the observed continuum emission may arise from 
multiple sources unresolved by the SMA beam; this issue is discussed further in §4.2. The 
total flux density of the three compact mm sources is 0.79±0.16 Jy; stated errors include the 
20% uncertainty in flux calibration. Thi s total corresponds to 15±3% of the single-dish flux 
density measured by IMinier et al.l (120051 ) with the SEST 15m telescope at 1.2mm (resolution 
24"). 

Figure [3] compares the morphology of the 1.3 mm dust continuum emission with the 
3.6 cm free-free continuum emission from the cometary UC Hll region, G192. 584-0. 041. 
Figure [3^1 shows the lower-resolution (l'/03 x 0"84) 3.6 cm VLA image superposed on the 
SMA 1.3 mm continuum, while Figure [3b shows the high-resolution VLA 3.6 cm image 
(0 / /27 x 0"23), with the positio ns of the newly-reported water m aser (see §3.2) and the Class 
I methanol masers detected by Kurtz. Hofner. fc Alvarez ( 2004 ) indicated. 



The ~ 1" resolution 3.6 cm VLA image presented in Figured provides the most detailed 
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view to date of the diffuse cometary "tail" of the UC H II region. The integrated flux density 
of G192. 584-0. 041 measured from this image is 26.0±0.1 mJy. Based on our measurement 
and published 2 cm integrated flux densities for G192. 584-0. 041 (IKurtz. Churchwell. fc Wood 



19941 ; iRengarajan fc Holll996l ). the spectral index from 3.6 cm to 2 cm is ~-0.1, consistent 



with optically thin free- free emission. The flux density is consistent with a single exciting 
star o f spectral type BO. 5 (as determined by lKurtz. Churchwell. fc Woodlll994l ; ISnell fc Bally 



19861 ). Extrapolating to 1.3 mm, we estimate the free-free contribution of G192. 584-0. 041 to 



the 1.3 mm flux density of SMAl to be <20 mJy (<3.5%). 

The 3.6 cm VLA image presented in Figure [3)d is the highest-resolution cm-wavelength 
image of G192. 584-0. 041 to date. With a resolution of 0'.'27 x 0'.'23, the continuum emission 
from G192. 584-0. 041 is resolved into three components (east to west): an arc, a point source, 
and an extended feature (the extended "tail" is resolved out in the higher resolution data). 
All three of these components overlap with the eastern side of SMAl, but none is coincident 
with the mm emission peak, in agreement with the estimate that the free-free contribution 
at 1.3 mm is quite small. The arc, which is oriented with its convex side towards the mm 
peak, contains the brightest 3.6 cm emission. The 3.6 cm peak is located in the southern 
part of the arc, east of the point source, and is offset by l'/l (~2800 AU) from the location 
of the SMAl mm continuum peak determined by fitting a single Gaussian component. The 
3.6 cm point source, which is located west of the arc and faces its concave side, is offset by l'/6 
(~4,200 AU) from the SMAl mm continuum peak. The peak brightness temperature of the 
3.6 cm point source is only 122 K at the current resolution. Absent a second radio frequency 
image with comparable resolution, it is not currently possible to ascertain the spectral indices 
of the individual components. The 3.6 cm images presented in Figure [3]^a-b) place strong 
limits on the presence of any additional Hn regions in S255N. Other than G192. 584-0. 041, 
no cm-wavelength emission is detected down to a 5<r limit of 90 /iJy beam -1 (high-resolution 
image) . 



3.2. Water maser emission 

Water maser emission was detected at the position 06 h 12 m 53 s .71, 18°00'27/6 (J2000), 
offset 0'/9 (~2,300 AU at 2.6 kpc) to the northeast of the 1.3 mm continuum emission peak 
of SMAl, as determined by fitting a single Gaussian component. This is the first report 
of water maser emission from S255N. The peak intensity is 2.8 Jy beam -1 (corrected for 
primary beam attenuation) at vlsr=+8.5 km s . The line is barely resolved by the 0.33 
km s -1 spectral resolution. 



The positions of the Class 1 44 GHz (7q-6i) A + methanol masers detected by 



Kurtz. Hofner. Alvarez 
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( 120041 ) are marked with crosses in Fig. [3b, wh ich shows the 1.3 mm and 3.6 cm continuum 
emission and the newly-reported water maser. iKurtz. Hofner. fc Alvarez! (120041 ) estimate an 
astrometric uncertainty of 0'/5 for the CH 3 OH maser spots, while the absolute astrometry 
of the H 2 maser is better than (X'l. The position of one of the CH 3 OH maser spots is 
consistent with SMAl, within the astrometric uncertainty. The newly detected water maser 
is < 1" from two CH3OH masers, and falls into the linear pattern of 44 GHz (7o-6i) A + 
CH3OH maser spots that extends northeast from SMAl. 



3.3. Line emission 

Molecular line emission from H 2 CO, CH 3 OH, SiO, CN, DCN, and HC 3 N is detected in 
S255N; the specific transitions, frequencies, and upper state energies are listed in Table [2j 
The lines detected i n S255N are the sam e as those detected in the spectral regions covered 



by our sidebands by ISutton et al.l (119851 ) in their line survey of the Orion A molecular cloud, 
which is similar to our data in spectral resolution (1.3 km/s), rms sensitivity (0.2 K), and 
linear size scale (30"= ~13,500 AU at 450 pc). Integrated intensity images for CH 3 OH, SiO, 
and H2CO are presented in Figure H(a-d) and for DCN, HC3N, and CN in Figure [5](a-c). 
The distributions of molecular emission observed in S255N fall into two main categories: 
CH3OH, H 2 CO and SiO exhibit emission from multiple locations, while DCN and HC3N are 
detected only in the vicinity of SMAl. CN exhibits compact emission towards SMAl, and 
is also weakly detected toward SMA3. The CN lines have the lowest E upper of the observed 
transitions, and the CN images show artifacts from large-scale emission resolved out by the 
interferometer, suggesting that much of the CN emission originates in an extended, cool 
envelope around the compact continuum sources. 

As shown in Figure IHa-d), the spatial distributions of the integrated emission from 
H2CO, CH3OH, and SiO are similar to one another. The kinematics of these molecules 
are complex, as illustrated in Figures [6] and Ui with multiple spatially and kinematically 
distinct components apparent. A finder chart for the positions of the profiles displayed in 
Fig. [7] (named after their relative positions with respect to SMAl) is shown in Figure [5(d). 
The positions are listed in Table [3j Line centroid velocities, Avfwhm, and integrated line 
intensities obtained from Gaussian fits to the line profiles at these positions are listed in Table 
HI Unless otherwise noted, the fit parameters in Table H] are for the strongest component in 
the spectrum. Fits to SiO line profiles are not included because the SiO line shapes are so 
complex. 

The strongest molecular emission in S255N lies toward the "SW" position ~ 6" to the 
southwest of SMAl at a peak velocity between ~ 6 — 8 km s _1 (Figs. IU El Et, Table HJ). 
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This molecular emission is not coincident with any mm continuum emission and overlaps 
the southern edge of the extended "tail" of the UC Hll region (Fig. H]). The line emission 
toward the "SW" position is broader than at any other location in S255N, with Avfwhm ~ 
7 and 9 km s _1 for CH 3 OH and H 2 CO, respectively, and shows pronounced blue wings. The 
velocity of the H2CO peak is slightly more blueshifted than CH3OH, while SiO is significantly 
blueshifted relative to both H 2 CO and CH 3 OH (Fig. [7b, Table 0}. 

The second-brightest region of H 2 CO, CH 3 OH, and SiO emission in S255N is located in 
the vicinity of SMAl, east of the cometary head of the UCHll region (Fig. HJ). DCN, CN, 
and HC3N have their strongest emission in this area (Fig. [5]). Spectral line profiles for DCN, 
CN, and HC3N are shown in Figure [8J and channel maps for DCN are shown in Figure [9j 
Two positionally and kinematically distinct components are evident in H 2 CO, DCN, CN, 
HC 3 N, and weakly, CH 3 OH, in the vicinity of SMAl. One component (denoted SMAl-NE) 
lies 1'.'21 to the northeast of the SMAl mm peak at a velocity of ~ 7 km s _1 , and the other 
(denoted SMAl-SW) lies T/17 southwest of the SMAl mm peak at ~ 11.5 km s" 1 (Figs. El 
[?k,b, and[8k,b). SMAl-SW is l'/08 east of the 3.6 cm point source. 

Interestingly, SMAl-NE and SMAl-SW also show differences in their chemical proper- 
ties. For example, H 2 CO shows nearly equal strength towards both positions, as does HC3N, 
while CH3OH is much stronger toward SMAl-SW (Fig. [6] [8]). In contrast, DCN and CN are 
both significantly stronger toward SMAl-NE (Fig. [8]). Some differences in the peak velocities 
at the two positions are also apparent amongst species. Relative to the other molecules, SiO 
is significantly blueshifted (v LSR < 5 km s" 1 ) towards both SMAl-SW and SMAl-NE (Fig.[7J 
Table HJ). DCN is slightly redshifted relative to CN and HC3N toward both positions (Fig. [HI 
Tabled]). The CN and HC3N lines are also more than twice as broad as those of H 2 CO or 
CH3OH toward SMAl-SW (Table ij. 

The H 2 CO, CH3OH, and SiO integrated intensity peak located ~5" north of the mm 
continuum source SMA2 (Fig. [U "NW" position in Fig. [5U), is comprised of relatively weak, 
broad emission (Fig. [7J). The H 2 CO and CH3OH lines are narrower than at the SW position, 
but broader than at any of the other positions (Tabled]). As at the other positions, the peak 
of the SiO line profile is blueshifted relative to H 2 CO and CH 3 OH (Fig. [7f). 

The line emission north and northeast of the SMAl region ("N" and "NE" positions, 
finding chart Fig. [5H) consists of narrow velocity features in CH3OH and H 2 CO, and, at 
the NE position, SiO (Figs. [6] & [7U-e). In contrast to the other positions, at the NE 
position CH3OH, H 2 CO, and SiO have the same velocity, vlsr~8 km s _1 (Fig. [7U, Ta- 
ble H]). The velocity of the H 2 CO and CH 3 OH emission at the N position is similar to 
the v LSR ~11.5 km s" 1 component toward SMAl-SW (Table ij). The H 2 CO and CH 3 OH 
lines are narrow, and the CH3OH peak is slightly blueshifted relative to H 2 CO. Broad and 
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weak blueshifted SiO emission is also detected at this position (Table HI Fig. [Tfe). 



3.4. Spitzer Space Telescope IRAC Observations 

Figure fTUl shows a three-color IRAC image (red 8.0 /im, green 4.5 /im, blue 3.6 /im) of 
S255N, overlaid with contours of the 1.3 mm continuum emission (yellow) and the 3.6 cm 
continuum emission (whi t e). The positions of the Class I CH 3 OH masers reported by 



Kurtz. Hofner. Sz Alvarea (120041 ) and the newly-reported water maser are marked with crosses. 
Most of the observed mid-IR emission is offset to the northwest of the UC H II region, and this 
diffuse emission appears in all IRAC bands. An exception is the linear, green 4.5 /im emis- 
sion feature that extends NE from the SMAl mm continuum peak. No mid-IR emission is 
associated with either SMA2 or SMA3, indeed these two positions are notably absent of IR 



emission. 



4. Discussion 



4.1. Mass estimates from the dust emission 



With an estimate of the dust temperature, we can estimate the masses of the compact 
dust sources S MAl, SMA2, and SMA3 using a simple isothermal model of optically thin 
dust emission (IBeltran et al.ll2006l ): 



M, 



gas,thin 



R F v D 2 
B(u,T d ) k v 



where R is the gas-to-dust mass ratio (assumed to be 100), F v is the observed flux density, 
D is the distance to the source, B(z/, T^) is the Planck function, and k u is the dust mass 
opacity coefficient. At 1.3 mm, the value of k for gas densities of 10 6 -10 8 cm -3 does not 
differ much for grains with thick or t hin ice mantles; we adopt a v alue of Ki.3 mm =l cm 2 g -1 
for all of the compact mm sources (lOssenkopf & Henning Ill994l ). The assumption of low 
optical depth is justified by the low observed millimeter brightness temperatures (Tabled]), 
however, for highest accuracy we have made the small correction to our derived masses for 
non-zero optical depth using the formula: M gas = M gaS)t hi n T/ (1 — e~ r ). 

Previous determinations of the dust temperature and mass in S255N have relied on 
fitting multiple components to the (unresolved) mid-IR to mm SED. iMinier et al.1 (120051 ) fit 
a hot, compact core (T = 106 K, diameter = 400 AU) and an extended warm envelope (T = 
44 K, diameter = 58,000AU) to a SED comprised of MSX, IRAS, SCUBA, and SEST data, 
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assuming a distance of 2.6 kpc. The derived luminosity and gas mass are 1.1 x 10 5 and 
220 M , respectively. While many of the datapoints in the SED constructed by lMinier et al. 



(120051 ) blend emission from all three SMA sources, the very compact hot core implied by 
their fits would be unresolved by our SMA beam (~12,200 x 6,200 AU at 2.6 kpc). Thus, the 
hot core temperature of 106 K derived from the SED modeling provides an upper limit for 
the dust temperature of the compact SMA sources. The fitted warm envelope temperature 
of 44 K is likewise a good lower limit to the temperature of SMAl since it dominates the 
1.3 mm flux, contributing 73% of the total SMA flux density. 

Several single dish estimates of the gas temperature are also available. For example, 
Effelsberg 100-m observations of NH3 (1,1) and (2,2) (resolution 40") suggest that the kineti c 
temperature of the gas is only T^ in = 23 ± 1 K (IZinchenko. Henning. fc Schreyerl 119971 ). 
Me asurements of CH, s C ? H (6-5) K=0-3 toward S255N with the Onsala 20-m (resolution 38") 
by iMalafeev et al.l (120051 ) yield T rot = 35 ± 1 K, in better agreement with the extended 
warm component derived for the dust. These authors find significantly higher temperatures 
using CH 3 C2H(6-5) than NH 3 towards all five sources observed (including S255) and suggest 
that methyl acetylene may preferentially trace warmer/denser gas. In any case, since beam 
dilution may play a significant role in these single dish estimates, they can only provide a 
lower limit to the gas temperatures on SMA sizescales. 

From our SMA line data it is clear that SMAl is the warmest of the compact mm 
sources: the two detected transitions with the highest E upper (both HC 3 N, E upper =131 K 
and 142 K) are detected only t owards SMAl, suggesting this may be a hot core (e.g. 
Hatchell. Millar, fc Rodgerslll998l ). DCN is also seen at this warm position. Though DCN is 
formed in cold clouds, in this case it can serve as a young hot core tracer since its presence in 
this warm region suggests it has recently been liberated from the icy mantles of dust grains 



^e.g. iMangum et al.lll99ll ). The HC3N emission is consistent with an upper temperature 
limit of ~100 K; a more quantitative determination is not possible with only two observed 
transitions. In contrast, the ratio of the H 2 CO(3o,3-2o,2) and H2CO(32,2-2 2 ,i) lines is a reli- 
able density-independe nt temperature diagnostic fo r Tk < 50 K, and A^(para-H 2 CO)/Af < 



10 13 - 5 cm 2 (km s 1 ) 1 ( IMangum fc Woottenlll993l ). In this regime, the 30,3-20,2/32,2-22,1 ra- 



tio ranges from 15 to 5 for Tfc=20 to 50 K. In contrast, the observed line ratios in S255N 
are less than 2.5 throughout the imaged region and are smallest (~ 1.5) toward SMAl, 
suggesting that the column density (i.e. opacity) and/or temperature is too high for these 
lines to be diagnostic. For the position of SMA1-NE, assuming a temperature of 75 K, 
N H2 co~3.7 x 10 13 cm" 2 from the H 2 CO(3o, 3 -2o,2) line and N H2 co~l-0 x 10 14 cm~ 2 from 
the H 2 CO(32,2-2 2 ,i). This comparison suggests that the 3 ,3-2 , 2 line is moderately optically 
thick compared to the 32,2~22,i line, and that the low line ratios are a combination of both 
the column density and temperature being higher than the diagnostic range of these two 
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transitions. Combining this analysis with the SED models and the single dish line results 
described above, the allowed temperature range for SMA1 is 40 - 100 K. The resulting ranges 
of gas mass, column density, and number density computed for SMAl are shown in Table 
The high derived gas density (n# 2 ~3-16xl0 6 cm -3 ), also implied by the presence of the 
water maser, indicates that the gas and dust temperatures are likely to be well- coupled (e.g. 



Kaufman. Hollenbach. fc Tielens 1 1 19981 : ICeccarelli. Hollenbach fc Tielens Ill996l ). 



Unlike SMAl, SMA2 and SMA3 are not accompanied by significant line emission. 
H 2 CO, CH 3 OH, and SiO emission are present to the north of SMA2, and CN emission 
is detected toward SMA3 (§3.3), but the physical relationship (if any) between this line 
emission and the dust continuum sources is unclear. Also unlike SMAl, SMA2 and SMA3 
are not associated with mid-IR emission in any IRAC band (§3.4). Instead, SMA2 and SMA3 
appear to be cold, dark, young mm cores, without evidence for current star formation. On 
the basis of the lack of line and mid-IR emission towards SMA2 and SMA3, we adopt a 
lower temperature limit of 20 K and an upper temperature limit of 40 K for these sources. 
The corresponding range of masses, column densities, and number densities for SMA2 and 
SMA3 are tabulated in Table El The mass of each (7-17 M for SMA2, 6-13 M for SMA3) 
is sufficient to form a low to intermediate mass star. No other cores are detected in the field 
to a 5a upper limit of M < 3M (at T = 20 K). 



4.2. Velocity Structure and Outflows 

Figure [10] shows a close-up view of the Spitzer 3-color IRAC image shown in Figure [TJ 
The brightest mid-IR emission is extended along a NE-SW axis, approximately parallel to 
the axis of the UCHll region, but with an offset to the northwest of ~ 2". The overall 
morphology of S255N is consistent with the multi-band bright mid-IR emission tracing the 
surface of the UCHll region, which is less dense to the southwest (of SMAl), as indicated 
by the diffuse "tail" of the cometary UC Hll region extending in this direction. However, 
the detailed interpretation of the mid-IR emission toward S255N is complicated by the offset 
described above, and the sharp cutoff of the 3.6 and 8.0 /zm emission along the southeast 
boundary of the UCHll region. Indeed, mid-IR emission is notably absent toward SMA2 
and SMA3, as well as toward much of SMAl. A likely scenario for this behavior is absorption 
of the mid-IR emission by the high column density mm cores; in this scenario the bulk of 
the relatively cold mm cores must be in front of the UCH II region. 

With the exception of SMAl, the molecular emission in S255N is not obviously associ- 
ated with any continuum emission, and is therefore unlikely to be centrally heated. However, 
as described in §4.1 the H 2 CO line ratios suggest the gas is warm. Thus, it is likely that 
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much of this emission is associated with outflow material, although the number of outflows 
and their driving source(s) are unclear. Published data on large-scale outflows in th e region 



|e.g. iMiralles et al.lll997l ; iRichardson et al.lll985l ; iHeyer et al.lll989l ; iRuiz et al.lll992l ) are un- 
fortunately too low in angular resolution to be useful in distinguishing outflows associated 
with S255N from those associated with S255IR to the south, and/or the maps are swamped 
by emission from a large outflow flowing north from S255IR. Excluding the "SW" position, 
the relatively narrow linewidths of these S255N line emission regions suggest that they are 
density enhancements within a larger extended flow resolved out by the interferometer. The 
relative similarity of the line center velocities further suggests that the outflows are mostly 
in the plane of the sky. 

The linear morphology of the (green) 4.5 /mi emission northeast of SMA1 is suggestive 
of an outflow (Fig. [TO]). Such 4.5 /mi nebulosity is a conspicuous feature of IRAC images 
of star forming regions. Recent analysis of the massive DR21 outflow, the best-studied ex- 
ample, has shown that H2 line emission accounts for ~50% of the observed 4.5 /mi IRAC 



(Davis et al. 


2007; 


Smith et al. 


2006) 



In S255N, the 

2.122 /mi H 2 clump S255:H2-3 lies at the base of the 4.5 /im nebulosity; the H 2 clump is 
also coincident, within reported astrometric uncertainties, with our newly-reported water 
maser. Both the 2 122 /im H 2 1-0 S(l) line and the 4.6947 /mi H 2 0-0 S(9) line, identified by 



Smith et al.l (|2006f) as the dominant contrib utor to IRAC band 2, trace moderate- veloc ity 
shocks (braine. Roberge. &; Dalgarn 

3E2S3). Recent models by ISmith fc Rose"nl (12005^ of 



shocks in dense protostellar molecular jets predict that the integrated H 2 line contribution 
to IRAC band 2 will be 5-14 times greater than to IRAC band 1 (3.6 /mi), consistent with 
the ratio of the emission seen in these bands toward the linear 4.5 /im feature. The 44 
GHz Class I methanol masers, five of which lie along t he 4.5 /mi emission featur e , pro - 



vide further evidence for its identification as an outflow. iKurtz. Hofner. fc Alvarez! (120041 ) 
found that masers of this type are often found in association with such outflow tracers 
as SiO (IRAS 20126+4104, G31.41+0.31 and G34.26+0.15) and H 2 (IRAS 20126+4104). 
In G31.41 +0.31, the 44 GHz methanol m asers ar e also associated wi t h ther mal methanol 
emission (IKurtz. Hofner. fc Alvarezl 120041 ) , which iLiechti fc Walmsleyl (119971 ) found traced 
shock/clump interfaces in the DR21 outflow. The parallels between these examples and 
S255N strongly suggest that the 44 GHz methanol masers, H 2 CO, CH3OH, SiO, H 2 , and 
4.5 /mi emission extending northeast from SMAl-NE trace a molecular outflow from a pro- 
tostar, probably SMAl-NE. 

At the SW position, the broad lines with strong blue wings combined with the morphol- 
ogy of the H 2 CO and CH 3 OH channel maps are consistent with a blueshifted outflow lobe 
driven by SMAl-SW. Notably, no 44 GHz methanol masers coincide with the very strong 
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thermal CH3OH emission of the SW line peak, although elsewhere in S255N the methanol 
maser flux densities are loosely correlated with the strength of thermal CH 3 OH emission. 
The absence of masers towards the SW position suggests that the phys ical conditions are 



not appropriate for the co llisional pumping of Class I methanol masers (ICragg et al.lll992 



Plambeck & Mentenlll990f ). 



4.3. The Nature of SMA1 

The complex kinematic behavior of the molecular line emission in the vicinity of SMAl 
including SMAl-NE, SMAl-SW, and position "SW" (§3.3) is difficult to explain in the 
context of a single protostar. Though SMAl-NE and SMAl-SW could be interpreted as the 
blue and red-shifted lobes, respectively, of an outflow, this scenario does not explain the very 
blue-shifted emission further to the southwest at position "SW" . Rotation also seems like an 
unlikely explanation for the velocity gradient between SMAl-NE and SMAl-SW since the 
gradient is parallel to the direction of the two probable outflow regions: the 4.5 fim emission 
to the northeast and "SW" to the southwest. Instead, the combination of the chemical and 
kinematic differentiation between SMAl-SW and SMAl-NE suggests the presence of two 
individual sources, one at +7 km s _1 and one at +11.5 km s -1 . To investigate this possibility, 
in Figure [TT] we show a uniform weighted SMA millimeter continuum image restored with 
a beam of l'.'O (~ 3 times smaller than the longest observed baseline), which essentially 
reveals the location of the clean components. The localization of clean components into two 
main regions in the vicinity of SMAl suggests the presence of at least two sources separated 
by ~l'/84 (4800 AU). If these two clean component peaks correspond to real dust sources, 
their positions are in good agreement with the two kinematically distinct formaldehyde peaks 
(< O'.'l and < 0'/5). The northeast component of the pair is also within 0'/2 of the water maser. 
That two distinct protostars would exist with this separa tion is reasonable, as multiplicity o f 



protostars has been observed on scales of <6,000 AU (e.g. lMegeath. Wilson, fc Corbinll2005l ). 
Although the presence of two protostars is a plausible interpretation, it clearly requires higher 
resolution continuum observations for confirmation. 

In any case, the molecular line emission from SMAl-NE and SMAl-SW is reminiscent of 
hot molecular cores (HMCs), particularly the detection of DCN and HC3N. These molecules- 
like CH3OH and n 2 CO, which also show strong emission towards SMAl-NE and SMAl-SW- 
are present in the gas phase in HMCs becau se they have been evaporated from grain mantles 
(e.g. ICasellill2005l ; ISzczepanski et al. 1120051 ). Complex organic molecules such as HCOOCH3, 



however, are believed to be "daughter" specie s, formed in the gas phase by reactions of 



"parent" species such as H 2 CO and CH 3 OH flCasellil 120051 ). Thus, while SMAl-NE and 
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SM Al-SW do not exhibit the truly cop i ous molecular emis sion observed towards some HMCs 



[e.g. lHatchell. Millar, fc Rodgerslll998l : ISchilke et al.ll2006[ ). this is consistent with SMAl-NE 
and SMAl-SW being very young sources, in which gas-phase hot-core chemistry has not yet 
produced abundant complex organic molecules. 

The line emission from the SMAl sources is unusual in that the DCN (3-2) emi s sion is 
stronger than HC3N(24-23). By modeling deuterium chemistry, iRoberts fc Millar I (120001 ) 
find that the steady-state abundance of DCN in molecular clouds is a complicated function 
of temperature and density (see their Figure 7), but generally higher at low metallicity. 
We note this because S255N is located approximately (I ~ 192°) in the direction of the 
Galactic anticenter, and may have lower meta llicity than inner-galaxy star-forming regions 
(jDaflon & Cunhal 120041 ; lAfflerbach et al.l 119971 ) . Our SMA data, however, do not allow us 
to disentangle the effects of abundance and excitation on the strengths of DCN and HC3N 



emission. 



The geometry of a HMC located a few arcseconds ahead of the vertex of a cometary 
UCHll region has been seen in other objects observed at high angular resolution, and this 
notable configuration has led to much discussion on the energy source resp onsible for these 
HMCs . Th e best-studied case for e xternal heating is G34.26+0.15, in which IWatt &: Mundy 
Jl999h and iMookerjea et al.l (120071 ) argue that HMC emission (characterized by complex ni- 
trogen and oxygen-rich molecules) arises in gas heated by component C, the most evolved 
of three nearby UCHll regions. In the absence of extinction, a 1.1 x 10 4 L Q source at 
the location of the 3.6 cm point source could heat SMAl-NE to ~37 K, and SMAl-SW to 
~51 K, consistent only with the lower end of the range of plausible gas temperatures. In 
contrast, G29.96-0.02 is the prototype for a HMC located ahead of a cometary UCH II region 



and internally heated by a high m ass protostellar object (IDe Buizer. Osorio. fc Calvetll2005 



Gibb. Wvrowski. & Mundvl 120031 ). In G29.96-0.02, HMC emission is coincident with a re- 
solved 1.4 mm continuum source, water maser spots, and a mid-IR sub- arcsecond point so urce 
(foe Buizer. Osorio. fc Calvetl I2OO5I : bibb. Wvrowski. fc Mundvl I2OO3I : loimi et al.l liooil . and 
references therein). In S255N, the arrangement of the 1.3 mm and 3.6 cm sources along 
with the presence of water maser emission coincident with the molecular emission closely 
resembles the case of G29.96-0.02. We favor the interpretation that one or more sources 
younger than the excitation source of the UCH II region are present and responsible for the 
compact dust and molecular line emission from SMAl, consistent with the interpretation of 
the hot core emission outlined above. 
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5. Conclusions 

Our multiwavelength observations of S255N reveal significant new details in this lu- 
minous star-forming region. While the previously-identified UCH II region dominates the 
cm continuum and mid-IR emission, the 1.3 mm continuum emission has been resolved into 
three compact cores (SMAl, SMA2, and SMA3) clustered on scales of 0.1-0.2 pc. Dominated 
by dust emission, these cores range in mass from 6 to 35 M Q . There are no mid-infrared 
point source counterparts to any of the dust cores, suggesting an early evolutionary phase. 
The spectral line emission at the position of the brightest core, SMAl, is spatially compact 
and includes HC 3 N, CN, DCN, CH 3 OH, SiO, and H 2 CO. SMAl appears to be a developing 
hot core offset by a few thousand AU from the UCH II region. The chemical and kinematic 
structure toward SMAl is suggestive of further multiplicity at these scales. A 4.5 /im linear 
feature emanating to the northeast of SMAl is aligned with a cluster of methanol masers and 
likely traces a outflow from a protostar within SMAl. We conclude that S255N is actively 
forming a cluster of intermediate to high-mass stars. In addition, we speculate that some of 
the missing flux in the SMA continuum image could be in the form of additional compact 
low-mass, cold dust cores that lie below the sensitivity limit of our observations (M ~ 3M 
at T = 20 K). Higher- resolution and more sensitive observations are needed to search for ad- 
ditional protostars in S255N and other young protoclusters. Resolving individual protostars 
in regions like these is a necessary task to determine how dense these young protoclusters 
are, and how interactions among protostars in protoclusters may affect the process of star 
formation. 

This work is based in part on observations made with the Spitzer Space Telescope, which 
is operated by the Jet Propulsion Laboratory, California Institute of Technology under a 
contract with NASA. This research has made use of NASA's Astrophysics Data System 
Bibliographic Services and the SIMBAD database operated at CDS, Strasbourg, France. 
Research at the CSO is funded by the NSF under contract AST96-15025. CJC is supported 
by a National Science Foundation Graduate Research Fellowship and acknowledges partial 
support from a Wisconsin Space Grant Graduate Fellowship. CJC would like to thank the 
SMA and NRAO for student research support. 
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Table 1. Properties of millimeter continuum sources in S255N 



Source 


J2000 coordinates 


A.3mm 


Size a 


rr b 

1 1.3mm 






a ( h m s ) 5 (° ' ") 


(Jy/b) 


[ " x " (•)] 


(Jy) 


(K) 


SMAl 
SMA2 
SMA3 


06 12 53.67 +18 00 26.9 
06 12 52.97 +18 00 31.9 
06 12 53.69 +18 00 18.5 


0.29 
0.09 
0.05 


3.9x2.0 (17.4) 
2.2x1.6 (97.5) 
5.2x <1.3 (169.6) 


0.58 + 0.12 
0.12 + 0.02 
0.09 + 0.02 


1.84 
0.82 
0.33 


Total 








0.79 + 0.16 





deconvolved source size determined by fitting a single Gaussian component to each 
source. The SMA beam is 4'/7 x 2'/4 (P.A.=-45.65°). 

b Uncertainties include 20% calibration uncertainty. 

c Brightness temperature computed using the Rayleigh- Jeans approximation. 



Table 2. Molecular species and transitions observed in S255N 



Species 


Transition 


Frequency E upper /k 






(GHz) (K) 



SiO 

DCN 

H 2 CO 

HC 3 N 

CH3OH 

H 2 CO 

H 2 CO 

CN a,6 
CN a,c 

CN a 
CN 
CN 
HC 3 N 



a These components are blended in our spectra. 

b The CN hyperfine components at frequencies 
lower than this one lie outside of our observed band- 
pass. 

c This transition is used to set the velocity scale 
for CN in Fig [8] 



5^4 
3^2 

3o,3 ~~ > 2o,2 

24^23 

4+2,2,0 ~~ *■ 3+1,2,0 



32,2 - ► 


2 2 ,i 


32,1 - > 


22,0 


2o,3,3 — 


> lo,2,2 


2o,3,4 — 


* lo,2,3 


2o,3,2 — 


* lo,2,l 


2o,3,2 — 


* 1q,2,2 


2o,3,3 — 


* lo,2,3 


25^24 





217.104980 


31.2 


217.238538 


20.9 


218.222192 


21.0 


218.324723 


131 


218.440050 


45.5 


218.475632 


68.1 


218.760066 


68.1 


226.874191 


16.4 


226.874781 


16.4 


226.875896 


16.4 


226.887420 


16.4 


226.892128 


16.4 


227.418905 


142 
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Table 3. Spectral line positions 



J2000 coordinates 



Name 


a ( h m s ) 


5 (° ' ") 


SMA1-NE 


06 12 53.73 


+ 18 00 27.8 


SMA1-SW 


06 12 53.64 


+18 00 25.4 


SW 


06 12 53.45 


+18 00 23.0 


NE 


06 12 53.76 


+18 00 33.2 


N 


06 12 53.70 


+18 00 39.8 


NW 


06 12 52.86 


+18 00 36.2 



Table 4. Fitted line properties 



H 2 CO(3q,3-2o,2) CH3OH DCN HC 3 N(24-23) 

Position Center Width / S v dv Center Width / S v dv Center Width / Svdv Center Width / S v dv 

km s — 1 km s -1 Jy b _1 *km s _1 km s -1 km s _1 Jy b _1 *km s _1 km s _1 km s — 1 Jy b _1 *km s~ 1 km s -1 km s _1 Jy b _1 *km s — 1 



SMA1-NE 


6.9(0.3) a 


5.9(0.6) a 


5.5(0.7) a 


6.6(0.3) c 


3.4(0.9) c 


0.9(0.3) c 


8.2(0.1) 


3.3(0.2) 


5.2(0.4) 


6.3(0.3) a 


3.5(0.8) a 


1.9(0.6)" 


SMA1-SW 


12.1(0.1) 


2.9(0.2) 


4.2(0.3) 


11.2(0.1) 


2.4(0.2) 


3.1(0.4) 


10.7(0. 3) b 


5.6(0.7) b 


3.1(0.5) b 


9.5(0.8) a 


11.2(2. l) a 


3.0(0.7) a 


SW 


6.3(0.1) 


9.2(0.4) 


24.3(1.2) 


7.8(0.1) 


6.9(0.3) 


13.4(0.8) 














NE 


8.3(0.1) 


2.9(0.3) 


2.4(0.4) 


8.3(0.2) 


2.5(0.4) 


1.6(0.4) 














N 


11.5(0.1) 


2.6(0.3) 


2.4(0.3) 


10.1(0.2) 


2.5(0.5) 


1.5(0.4) 














NW 


8.9(0.2) 


5.8(0.6) 


5.3(0.7) 


8.3(0.2) a 


5.7(0.5) a 


3.7(0.5) a 















a Not well fit by a single Gaussian. 

b Gaussian fit encompasses two blended components. 

c Paramctcrs for second-strongest velocity component. Strongest component is very similar to that towards SMA1-SW. 
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Table 5. Range of estimated masses of dust cores in S255N 





K 


T 

± dust 


7~dust 


M 


Nn 2 b 




Source 3, 


(cm 2 g -1 ) 


(K) 


(1.3mm) 


(M ) 


(10 23 cm- 2 ) 


(10 6 cm- 3 ) 


SMA1 C 


1 


40-100 


0.04-0.02 


35-13 


10.5-3.8 


15.9-5.8 


SMA2 


1 


20-40 


0.02-0.01 


17-7 


5.0-2.2 


7.6-3.3 


SMA3 


1 


20-40 


0.01-0.01 


13-6 


3.9-1.7 


5.9-2.5 


Total 








65-26 







a Assumed distance is 2.6 kpc. 

b Beam-averaged quantities.The SMA beam is 4'.'7 x 2'.'4 (P.A.=-45.65°). 

c The mass of SMAl was calculated using the 1.3 mm flux density less 
the estimated free-free contribution of 20 mJy. 
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Fig. 1. — Three-color Spitzer IRAC image of S255N and its surroundings showing 8.0 
fim (red), 4.5 /im (green), and 3.6 /mi (blue). S255N lies in a complex region of past 
and ongoing massive star-formation. 
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Fig. 2. — Greyscale and solid contours of the SMA 1.3 mm continuum with dotted contours 
of the CSO 350 jum continuum superposed. The primary beam of the SMA (56" at 218.7 
GHz) is indicated with a black circle. Naming conventions for mm and submm sources used 
in the literature and in this paper are also indicated. The black 1.3 mm contour levels 
are (-3, 3, 7, 15, 31, 47, 63) x 4 mJy beam -1 (the rms noise), observed with a 4"7 x 2"4 
(P.A.=-45.65°) beam. The dotted 350 fun contour levels are (2, 2.5, 3, 4, 5) x 16.5 Jy 
beam -1 , resolution 15". 
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Fig. 3. — (a) Black contours of the SMA 1.3 mm continuum, observed with a 4'/7 x 2'/4 
(P.A.=— 45.65°) beam, with greyscale and dotted contours of the 3.6 cm continuum, observed 
with an l'/03 x (X'84 beam (P. A. =—81.42°), superposed. The black 1.3 mm contour levels are 
(-3, 3, 7, 15, 31, 47, 63) x 4 mJy beam -1 (the rms noise). The dotted 3.6 cm contour levels 
are (-3, 3, 5, 7, 11, 21, 41, 61, 81, 101, 121, 161) x 34 /iJy beam -1 (the rms noise). The 
SMA beam (black ellipse) and VLA beam (filled black ellipse) are plotted at lower left, (b) 
Black contours of the SMA 1.3 mm continuum with greyscale of high-resolution (0'/27 x 0'/23, 
P.A. =77.51° beam) 3.6 cm continuum superposed. The positions of methanol masers and of 
the newly-reported water maser are marked. The VLA beam (filled black ellipse) is plotted 
at lower left. 
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Fig. 4. — (a-d) Colorscale of the 1.3 mm continuum with black integrated intensity contours 
for molecules that exhibit emission from multiple positions. The molecular species and 
upper state energies are indicated in the lower right of each panel (also see Table 2). The 
white contours show 3.6 cm emission from the UC H II region G192. 584-0. 041. The blue cross 
marks the position of the newly detected water maser. The black integrated intensity contour 
levels are (-3, 3, 5, 7, 9, 11, 13, 15, 17, 19, 21, 23, 25, 27, 29) x the rms noise levels: (a) 
CH 3 OH 0.82 Jy beam-^km s" 1 , (b) SiO 1.35 Jy beam'^km s _1 , (c) H 2 CO 3 , 3 -> 2 , 2 0.8 
Jy beam _1 *km s _1 , (d) H 2 CO 3 2;2 - * 2 2 ,i 0.76 Jy beam _1 *km s -1 . The white 3.6 cm contour 
levels are (-3, 3, 5, 7, 11, 21, 41, 61, 81, 101, 121, 161) x 34 /iJy beam^ 1 (the rms noise). The 
4'!7 x 2'.'4 (P.A.=-45.65°) SMA beam (filled black ellipse) and 1'.'03 x 0'.'84 (P.A.=-81.42°) 
VLA beam (filled white ellipse) are shown at lower left in each panel. 
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Fig. 5. — (a-c) Similar to Figure H^-d except integrated intensity images are shown for 
molecules that are only detected in the vicinity of SMAl. The black integrated intensity 
contour levels are (-3, 3, 5, 7, 9, 11) x the rms noise levels: (a) CN 0.9 Jy beam _1 *km s _1 , 
(b) DCN 0.55 Jy beam _1 *kin s _1 , and (c) HC 3 N 24^23 0.68 Jy beam _1 *kni s -1 . (d) Finding 
chart for representative line profiles shown in Figures [7] & [8] and discussed in §3.3. 
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Fig. 6. — Channel maps of (a) H 2 CO(3o,3-2 0i2 ) and (b) CH 3 0H, showing line emission 
(contours) overlaid on the 1.3 mm continuum (greyscale). The contours levels are (1, 2, 3, 
4, 5, 6, 7, 8, 9) x 0.28 Jy beam -1 . Each panel is labeled with the channel velocity. The 
position of the 3.6 cm point source is marked with a white cross. 
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Fig. 7. — Representative line profiles for molecules that exhibit emission from multiple 
positions, demonstrating the wide range of velocity and chemical behavior observed. The 
positions for which spectra are presented are indicated on Figure EH and listed in Table El 
A vertical line is drawn at 10 km s -1 for reference. 
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V |sr (km s V |sr (km s 1 ) 

Fig. 8. — Representative line profiles for molecules that have their strongest emission towards 
SMAl. The positions for which spectra are presented are indicated on Figure EH and listed 
in Table El A vertical line is drawn at 10 km s™ 1 for reference. Note that the vertical scale 
is not the same as in Figure [71 For CN, the weaker features offset by -16 and -22 km s _1 
from the main feature are due to the hyperfine components (see Table [2]). 
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Fig. 9. — Channel maps of DCN emission (contours) overlaid on the 1.3 mm continuum 
(greyscale). Each panel is labeled with the channel velocity. The position of the 3.6 cm 
point source is marked with a white cross. 
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Fig. 10. — Close-up three-color Spitzer IRAC image of S255N showing mid-IR emission 
offset to the NW of the UC Hn region with yellow SMA 1.3 mm and black 3.6 cm contours 
superposed. The colorscale correspond to: 8.0 /zm (red), 4.5 jttm (green), and 3.6 jttm (blue). 
The yellow SMA 1.3 mm continuum contour levels are (3, 7, 15, 31, 47, 63) x 4 mJy beam -1 
(the rms noise). The black VLA 3.6 cm continuum contour levels are (3, 5, 7, 11, 21, 41, 61, 
81, 1 01, 121, 161) x 34 /iJy b e am -1 (the rms noise). Positions of Class I CH 3 OH masers 
from iKurtz. Hofner. fc Alvarea (120041 ) are marked with red crosses, and the position of the 
newly-reported water maser is marked with a blue cross. 
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Fig. 11. — Greyscale image of the high resolution VLA 3.6 cm emission shown in Figure [3b 
with SMA 1.3 mm uniform weighted continuum contours superposed. The 1.3 mm image 
was restored with a 1" beam (~ 3 times smaller than the longest observed baseline) which 
essentially shows regions where clean components are concentrated. Red + symbols show 
the 44 GHz methanol maser positions, the blue A shows the H2O maser position, the green 
x symbols show the peak locations of the ~ 7 (SMA1-NE) and ~ 11.5 km s _1 (SMAl-SW) 
H 2 CO components, and the black o shows the location of SMA1 reported in Table [TJ The 
localization of clean components into two distinct regions in the vicinity of SMA1 suggests 
the presence of at least two continuum sources; this result requires confirmation by higher 
resolution mm/submm data. 



